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Abstract 

We propose the phase dynamics induced by spin waves in a superconduc- 
tor/ferromagnet /superconductor (SC/FM/SC) Josephson junction. The resistively shunted 
junction (RSJ) model, which describes the dynamics of superconducting phase difference, is 
extended to include the spin wave excitation by ferromagnetic resonance (FMR) using the 
gauge invariant phase difference between two s-wave superconductors. The current-voltage 
characteristics show step structures when the magnetization in FM is driven by tuning the 
microwave frequency to FMR in the SC/FM/SC junction. The result presents a new route to 
observe the spin wave excitation using the Josephson effect. 

PACS numbers: 74.50. +r, 76.50.+g 
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The Josephson effect is characterized by a zero voltage current through a thin insu- 
lating barrier sandwiched by two superconductors [1]. This effect is a macroscopic quan- 
tum phenomenon involving phase coherence between two superconductors. When a finite 
voltage-drop (V) appears in the junction, the phase difference (8) evolves according to 
86/ dt = 2eV/h. As a result, an alternating current flows due to the time dependence of 8. 
Moreover, when an alternating electric field due to microwave irradiation is applied to the 
junction, the current- voltage (I-V) characteristics show Shapiro stepsj^L Similar phenomena 
are observed in a normal metal weak link due to the proximity effect 3|, |4| . The Josephson 
junctions are useful devices. For instance, the ac Josephson current is utilized to detect the 
electron-spin resonance in the Josephson junction doped with magnetic impurities [5]. The 
dynamics of the Josephson junction which are reflected in the I-V characteristics are de- 
scribed by the resistively shunted junction (RSJ) model, which is used to analyze Josephson 
dynamics 0, Q|. 

The Josephson effect in superconductor/ ferromagnet /superconductor (SC/FM/SC) junc- 



tions has been of considerable interest in recent years[8|, [9|, |10|, lUl, ll2|, ll3|, |14|, 115J. One of 
the most interesting phenomena is the oscillation of the Josephson critical current (I c ) in 
SC/FM/SC junctions. The origin of the oscillation is the exchange splitting of the conduc- 
tion band in the FM, in which Cooper pairs penetrate into the FM with a finite center of 
mass momentum. If the thickness of the FM is about half of the period of the oscillation, 
the current-phase relation is shifted by rr from that of a conventional Josephson junction 
(0-iunction). This is called a 7r-junction, which has a potential application as a quantum 
bit 16j. The 7r-junction is also used to measure experimentally a nonsinusoidal current-phase 



relation 



171 ]. Most of the studies on SC/FM/SC junctions have so far been focused on the 



effect of the exchange splitting on Cooper pairs penetrating into the FM. Ferromagnetic 
materials possess dynamic properties such as spin waves, which may be excited by using 
ferromagnetic resonance (FMR). Therefore, the interaction between Cooper pairs and spin 
waves in the FM is expected to play an important role in the transport. Up to now, the 
spin dynamics in the FM has been disregarded in SC/FM/SC junctions. 

In this Letter, we study the effect of the spin dynamics on the Josephson effect in a 
SC/FM/SC junction. The RSJ model is extended to include the effect of spin wave excita- 
tions using the gauge invariant phase difference between superconducting leads. We adopt a 
model in which the magnetization of FM exhibits a precessional motion. The I-V character- 
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istic is calculated using the extended RSJ model. It is found that this characteristic exhibits 
step structures and the dc Josephson current is induced due to the spin wave excitation by 
FMR. The step structures appear whenever the relation uj = 2£Q is satisfied with i being 
an integer, where Q is the FMR frequency, and cuj = 2eV/h. We propose a new route to 
observe spin waves using the Josephson current in SC/FM/SC junctions. 

The system considered is a Josephson junction with a ferromagnet sandwiched by two 
s-wave superconductors. We choose a coordinate system such that the electrode surfaces 
are parallel to the yz and xz planes as shown in Fig. [TJ and the current flows along x. The 
Josephson junction is characterized by the phase difference 9(y,z,t) between SCs, which is 
described by the RSJ model as follows 

1 = *cSm^,,,t) + --^^ + C- - , (1) 

where i is an external current density and i° c is the critical current density. The flux quantum 
is denoted by $ - The resistance, R, and the capacitance, C, in the Josephson junction are 
normalized by the junction area, S yz , as R = R /S yz , C = C /S yz . R and C are the 
resistance and the capacitance in the Josephson junction, respectively. 

We consider the situation, in which the FM is exposed to a circularly polarized microwave. 
The microwave causes the precessional motion of the magnetization, which corresponds to 
the excitation of the uniform mode of a spin wave. The motion of magnetization due to the 



microwave radiation is described by the Landau-Lifshitz-Gilbert (LLG) equation 191 ] 



dM „ „ a 

-jM xH cS + — 



dM 
M x 

dt 



(2) 



dt M 

where M (t) is the magnetization of the FM, 7 is the gyromagnetic ratio, and a is the Gilbert 
damping. The effective field, to which M (t) responds, is given by if e g = H + fa ac (t), 
where H is an uniaxial magnetic anisotropic field, which is parallel to z axis and h ac = 
(/lac cos fit, h ac sin fit, 0) is the microwave driving field with frequency Q. In Eq. (j2J), we 
neglect an anisotropic precession of the magnetization due to the demagnetization field for 



simplicity 2Cj, |21( . The linearized solution of Eq. (j2j) is given by 



which describes the precession of M (t) around H with frequency Q and has a resonance 
with Q = jH , where M ± {t) = M x (t) ± iM y (t). The rotating magnetic field is induced in 



the FM through the magnetic flux density B (t) = AttM (t) . Due to the magnetic field, the 
superconducting phase is not gauge- invariant. Therefore, we derive the equations for the 
gauge invariant phase difference between SCs. In the SC/FM/SC junction, the magnetic 
field in the FM is equal to 4ttM (t), since we assume that the two superconductors separated 
by the FM of the thickness d are thick compared with the London's penetration depth and 
the magnetic field produced by the Josephson current is negligible. Therefore, the flux A$ y 
parallel to the y axis enclosed by a contour C in Fig. [1] is A$ y = 4irM y (t) d/^oAz. In a 
similar way, the flux A$ z parallel to the z axis is = 4nMd/^ Ay. Combining these 
relations with a gauge invariant phase difference, we obtain the differential equations with 



respect to 9 (y, z, t)|22|, |23(, 



4:7ld 

V y}Z e(y,z,t) = -—M(t)xn, (4) 

<4> 

where \J y , z = (0,d/dy,d/dz), n is the unit vector perpendicular to the yz plane. Equation 
(jlj) is satisfied by the following solution, 

9 (y, z,t) = 9 (t) - -—^y + ^-z, (5) 

where we adopt the gauge, such that the superconducting phase difference is equal to 9 (t) 
without magnetic field. Introducing Eq. (jSJ) in Eq. ([1]) and integrating over the junction 
area, we obtain the extended RSJ model that includes the spin dynamics in the FM, 



sin (*.M„ (t)J i0(T) d 2 0(T) 



with 



L z dAn 2 -fM z h ac 

S " $o Oo ' { ) 

j£ _ (n/n - 1) sin (^tjt) an/n cos (^tjt) 

v (n/n - if + ( a n/n ) 2 (o/o - 1) 2 + («n/n ) 2 ' 

where / is an external current and I c = I® (ir^o/^z) sin (7r<& 2 / $o). ^? ^ s the Josephson critical 
current, and $ z is the total magnetic flux parallel to the z-axis in the FM. The other variables 
are given by, r = t/Vj, tj = $ / (2nI c R), and /3 C = RC/tj. When $ z = and $ s M y (r) = 0, 
Eq. (jHJ) reduces to a conventional RSJ model. The extended RSJ model is applicable to both 
0- and 7r-junctions. Below, we solve numerically Eq. (J6j) for f3 c = (overdamped junction) 
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and calculate the I-V characteristic using the relation, (89/ dt) = 2e{V) /K where (• • •) 
denotes a time average. For the numerical calculation, a is taken to be 0.01 [l9|. 

First, we consider the situation that the magnetic field component of microwave is applied 
to the junction and the FMR occurs in the FM. This condition is realized by using a 
microwave cavity 2J]. The I-V characteristic is calculated using Eq. ([2D and the result 



is shown in the solid and dashed lines in Fig. 12(a). The vertical axis is the normalized 
current, I / I c , and the horizontal axis is the normalized voltage, V/I C R. As shown by the 
dashed line in Fig. [2(a), when the microwave frequency (O) deviates from the ferromagnetic 
resonance frequency (Oo), the I-V characteristic is in agreement with that of a conventional 
Josephson junction. On the other hand, when the microwave frequency is in the condition 
| 1 — Q/Q |< a at which FMR occurs in the FM, the I-V characteristic shows the step 
structure under the FMR as shown by the solid line in Fig. [2(a). 

Second, we consider the case that the electromagnetic field due to microwave irradiation 
is applied to the junction. For the electric field component, we use the ac current source 
model in which I is replaced by I + J ac sin (fit) in Eq. (jHJ), where 7 ac is the amplitude of ac 
current 25j. Figure [2(b) shows the I-V characteristic which reflects the electromagnetic field 
irradiation. In the dashed line in Fig. [2(b), when Q deviates from Q , the I-V characteristic 
shows Shapiro steps due to the ac current similarly to conventional Josephson junctions. 
When f2 is equal to that of the alternating Josephson current, these steps arise at V = n±Ml 
with n being an integer. When the microwave frequency is in the condition | 1 — Q/Qq |< 
a, the step structure dramatically changes as shown by the solid line in Fig. [2(b). The 
amplitude of the second step is much larger than that of the first step. This behavior is 
different from that of the Shapiro step for which the amplitude of the steps decreases with 
increasing the voltage. In addition, the step structures appear in the higher voltage region 
(V/I C R > 1.5) unlike the case of the Shapiro steps. 

To elucidate the origin of the step structure in the I-V characteristic due to the spin 
wave excitations, we analyze the first term of the Josephson current ij (t) in Eq. ([6]). Using 
the generating function of Bessel functions and the standard mathematical expansion of the 
sine function, we have 
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oo oo 

i[6(t)+2(n+m+l)Slt+2(n+m+l)ip] 
Cs " 



n=0 m=— oo 

oo oo 



+ ]T j 2m+1 ( Cs ) e -i[0(t)-2(n+m + l)m+ 2 (n+m+l)V] ) (g) 

Cs 



n=U m=— oo 



where J m (Cs) is the Bessel function of the first kind, Cs = &s/(Q — Qo), an d i> = 
arctan ( n °f^ o j ■ In Eq. ([9]), when 9 (t) = ±2 (n + m+1) Qt, the time averaged ij (t) has 
a nonzero component of the dc Josephson current. As a result, the I-V characteristics rep- 
resent the step structure without changing voltage. Using the relation, 86 (t) jdt = uj, we 
obtain the voltage at which the step structures occur, 

uj = 2{n + m+l)Q, (10) 

where n and m are integers. Figure [3] shows the amplitude of the step structure as a function 
of the frequency 2 (n + m + 1) Q in the I-V characteristic. The vertical axis is the normalized 
amplitude of the step. It is seen that the step structures only appear at even number unlike 
the case of the Shapiro step. 

Finally, we discuss the condition 2(n + m + 1)Q occuring step structures in the I-V 
characteristic. We convert Eq. (TTUT) into V = 2(n + m + l)^hQ, which is quite similar to the 
condition V = n^-Hft of the Shapiro step, where the voltage is proportional to the number 



of photon n and its energy Mlj^, |25|. In an analogous way, we interpret the above results 
as follows. In the present situation, a uniform mode of spin wave is excited in the FM by 
the FMR due to the microwave irradiation. Therefore, 2 (n + m + 1) and Ml in Eq. f fTUj) 
correspond to the number of quantized spin waves (magnons) and their energy, respectively. 
If Cooper pairs penetrating into the FM are scattered by the odd number of magnons, 
the Josephson coupling vanishes between the s-wave SCs because of the formation of a 



spin triplet state [26|, resulting in no Josephson current. On the other hand, the Josephson 
coupling occurs when Cooper pairs penetrating into the FM compose the spin singlet state 
because of the even number of magnon scattering, so that the Josephson current flows. 
Therefore, the even number of magnons couples to the Josephson current. 

In conclusion, we have studied the effect of the spin dynamics on the ac Josephson effect 
in the superconductor /ferromagnet /superconductor (SC/FM/SC) junction. The resistively 
shunted junction (RSJ) model is extended to include the spin wave excitations induced 



by the ferromagnetic resonance (FMR) and the equations for the gauge invariant phase 
difference between SCs have been derived. We find that the current-voltage characteristics 
show step structures when the relation uj = 2£Q with £ being an integer is satisfied, where 
ujj = leV I h and Q the FMR frequency. The results provide a new route to observe the spin 
waves using the Josephson current in SC/FM/SC junctions. 
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FIG. 1: Josephson junction with FM between two s-wave SCs. A microwave radiation /i ac onto 
FM causes the precession of the magnetization of FM. Coordinate system is chosen such that the 
electrodes surfaces are parallel to yz-plane. C is a contour in xz-plane. 
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FIG. 2: Current-voltage characteristics. The solid and dashed lines show the I-V characteristics 
microwave frequency for = 0.9 and Q/Qq = 1, respectively. We choose <£ s = 0.03. (a) 

Wlc = 0. (b) Jac/Je = 0.6. / ac is ac current across the junction due to the electromagnetic field. 
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FIG. 3: The amplitude of the step in the I-V characteristics vs the normalized Josephson fre- 
quency. The positions of step lie on the even number. Al2i/I c is the amplitude of step in the I-V 
characteristics, and I is an integer. 
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